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f^ ' Abstract. We rigorously prove a form of disorder-resistance for a class of one-dimensional 

^P . cellular automaton rules, including some that arise as boundary dynamics of two-dimensional 

solidification rules. Specifically, when started from a random initial seed on an interval of length 
L, with probability tending to one as L — >■ oo, the evolution is a replicator. That is, a region of 



o 

^r^ . space-time of density one is filled with a spatially and temporally periodic pattern, punctuated 

by a finite set of other finite patterns repeated at a fractal set of locations. On the other hand, 
the same rules exhibit provably more complex evolution from some seeds, while from other seeds 
their behavior is apparently chaotic. A principal tool is a new variant of percolation theory, in 
l2 ', the context of additive cellular automata from random initial states. 
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1 Introduction 

Cellular automata (CA) started from seeds, i.e., finite perturbations of a quiescent state, have 
been the subject of much empirical analysis, starting with [WoUj . The observed behavior falls 
roughly into four categories: (a) the perturbation remains localized in the sense that it never 
affects sites outside a bounded interval; (b) a periodic structure develops and spreads; (c) a 
replicating (also called nested or fractal) structure develops, with a recursive (but sometimes 
complicated) description; (d) unpredictable chaotic (or complex) growth generates a space-time 
configuration with apparent characteristics of random fields. Many CA are capable of behavior 
in multiple categories depending on the choice of seed, and this is true even for some of the very 
simplest one-dimensional CA. An example is the Exactly 1 rule, in which a cell is alive whenever 
exactly one of itself and its two neighbors were alive at the previous generation. Exactly 1 is 
capable of periodic, replicating, and chaotic behavior for different seeds; see |GG3] . 

If a particular CA is capable of chaotic behavior from some initial seed, it appears natural 
to conclude, by analogy with the second law of thermodynamics, that such behavior should be 
generic for that CA, in the sense that almost all sufficiently long seeds yield chaotic evolution. 
Shadowing results from dynamical systems [Pilj, with their general message of stability of chaotic 
trajectories, would also tend to support such a conclusion. Indeed, strong empirical evidence 
confirms that chaotic behavior is prevalent for many CA including Exactly 1; see |GG3] . 

In this article we exhibit a class of one-dimensional CA rules for which we rigorously prove 
that the opposite conclusion holds. Typical (random) long seeds self-organize into replicating 
structures, while exceptional seeds yield more complex behavior, including apparently chaotic 
evolution. 

We focus on one-dimensional range-2 CA rules with 3 states (although our techniques in 
principle apply to more general one-dimensional rules). Thus, the configuration of the CA 
at time t G {0,1,2,...} is an element ^t = {^t{x))xez of {0,1,2}^, and for a given initial 
configuration ^q, the evolution is given by 

^t+i{x) = f(Ct{x - 2),^t{x - l),6(x),6(x + l),Ct{x + 2 

for all x, t and a fixed function / (the CA rule). (In many cases the dependence on ^t will actually 
be restricted to the range-1 neighborhood x — 1, 2;,2;-|- 1). We sometimes write ^(x,t) = ^t{x) for 
the state of S, at the space-time point {x,t) € Z x [0,00). In keeping with standard convention, 
diagrams of space-time evolution are drawn with the space coordinate x increasing from left to 
right, and the time coordinate t increasing from top to bottom. 

A key supporting role will be played by the 1 Or 3 CA, a simple 2-state rule denoted by Aj, 
and defined as follows. The states are and 1, and the evolution is 

Xt+i{x) = Xt{x - 1) + Xtix) + Xt{x + 1) mod 2. 

As is well known, the additive structure of this rule enables many of its characteristics to be 
fully understood. (See Figure [TTT] below for an illustration.) 



1 INTRODUCTION 




Figure 1.1: Left: the configuration A* of 1 Or 3, started from a single occupied cell, up to time 
t = 32. Right: schematic depiction of a rephcator. The striped regions are fihed with a doubly 
periodic ether. The thickness of the white "buffer zones" remains constant for all time. 

We consider 3-state CA rules with the following special property. For any configuration ^, if 
we define Xt{x) = l[^t{x) = 1] for all x,t, then A evolves precisely according to the 1 Or 3 CA. 
We also assume that state is quiescent - that is, if .^o = then ^i = 0. We call any CA rule 
satisfying these two conditions a web CA. The idea is that the Is form an additive "web" which 
is not influenced by the distinction between Os and 2s, while the web may affect the 2s. As we 
will see later, web CA also arise in analysis of two-dimensional solidification CA. We will usually 
be interested in evolution from a seed, i.e. an initial configuration ^o with finite support. 

One of the simplest web CA, which we call Web-Xor, is defined by setting £,tix) = 2 if and 
only if At(x) = and there is a exactly one 2 among ^t-iix — l),S,t-iix + 1). (Together with 
the web CA condition, this is sufficient to specify the rule). Thus, 2s perform a 2-neighbor 
exclusive-or rule on the points that are not occupied by Is. Figure [L2] illustrates the evolution 
of Web-Xor from four different seeds. (States are always colored as: white; 1 black or grey; 
2 another color depending on the rule.) Our results imply that typical seeds result in behavior 
similar to the first picture. More specifically, we will prove that for certain classes of web CA, 
evolution from long random seeds yields with high probability a space-time configuration that is 
periodic except within some finite distance of an additive web. To state this conclusion precisely 
we need some more notation. 

An ether is an element rj of {0, 2}^ that is periodic in both coordinates. Two ethers 
are equivalent if one can be obtained from the other via some translation of Z^. In a CA 
configuration ^ we say that a set i^ C Z x [0, oo) is filled with r/ if ^ agrees with some ether 
equivalent to 77 on X. Let A* be the 1 Or 3 CA started from the seed consisting of a single 1 at 
the origin, and let A = {{x,t) : X*{x,t) = 1} be its support. See Figure [TTTl Let A(r) C Z^ be 
the set of space-time points at ^""^-distance at most r from A. 

For a given CA, we say that a seed ^0 (oi' equivalently the resulting configuration ^) is a 
replicator of thickness r and ether rj if each bounded component of Z^ \ A(r) is filled with rj. 
See Figure [LTJ It is a straightforward fact that A(r) has density as a subset of I? for any r. 
Therefore, in a replicator, the density of 2s within the cone {(x, t) : |x| < t} equals the density of 
the ether. Furthermore, it may be shown that for any replicator (of any CA), the configuration 
^ can be fully described in terms of a finite set of local patterns that are repeated at infinitely 
many locations prescribed by A. (This is the reason for the name replicator.) For more details 
we refer the reader to |GGP| . where the concept was introduced. 
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Figure 1.2: Four configurations of Web-Xor. The first (top) example, a replicator with zero 
ether, starts from a random string of 64 Is and 2s. The second and third examples, with 
respective seeds 12 and 11111012, are quasireplicators. The bottom example, with seed 1100112, 
is apparently chaotic. 
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Our results will apply to web CA rules satisfying two conditions which we call diagonal- 
compliance and wide- compliance. The conditions state that flow of information concerning the 
distinction between Os and 2s is blocked by certain local patterns of Is. The formal statements 
of the conditions are straightforward but somewhat technical, and we therefore postpone them 
to the next section. For now we note that Web-Xor is diagonal-compliant. 

A uniformly random binary seed on [0, L] is an initial configuration ^o in which (,o{x) 
takes values 0, 1 with equal probabilities independently for all x € [0, L], and outside [0, L]. 

Theorem 1.1 (Replication from random seeds). Consider a web CA that is either diagonal- 
compliant or wide- compliant, started from a uniformly random binary seed on [0,L]. There 
exist a random variable Rl taking values in [0, oo], and a random ether rjL (both deterministic 
functions of the seed), with the following properties. We have F{Rl = oo) — )• as L ^ oo, and 
indeed the sequence {Rl)l>o ^s tight. On the event Rl < oo, the configuration S^ is a replicator 
of thickness Rl + L and ether rjL. Furthermore, if any finite set of Os in ^q are changed into 2s, 
the same statement holds with the same Rl and r]L. 

Web CA rules may be further classified in the following way, which has implications for their 
production of ethers. A CA has no spontaneous birth (of 2s) if whenever ,^o contains no 
2s, .^1 also contains no 2s. Web-Xor has no spontaneous birth. Figure 11.31 shows four possible 
evolutions of a CA rule called Piggyback (to be defined in the next section) that is wide-compliant 
and has spontaneous birth. 

Theorem 1.2 (Trivial and non-trivial ethers). Assuming the conditions of Theorem \1.1[ Rl 
can be chosen to have the following additional properties. 

(i) // the CA rule has no spontaneous birth, then r]L = whenever Rl < oo. 
(ii) Suppose that the CA rule has spontaneous birth. If for some deterministic ether r] we have 

Rl < oo and tjl = r] for some binary seed, then for uniformly random binary seeds we 

have 

liminf P{Rl < oo and r]L = i]) > 0. 

Given any particular seed, there is a simple procedure to compute the random variable Rl 
appearing in Theorems 11.11 and 11.21 ai^d in particular to determine whether it is finite. (See 
Sections [8] and [9] for details). For many CA of interest, including Piggyback, there are multiple 
non-equivalent ethers r] for which the condition of Theorem ll.2( ii) indeed holds, and which hence 
have asymptotically non-trivial probabilities. The first two pictures in Figure 11.31 show two 
examples. Our methods allow the computation of explicit rigorous lower bounds on asymptotic 
probabilities of particular ethers. For example, in Piggyback, for the ether that results from the 
periodic initial state (00022222)°°, the liminf in the theorem is at least 0.1297, while (0)°°, (2)~ 
and (00002022)°° have lower bounds 0.5, 0.0398 and 0.0151 respectively. (In fact, more than 
100 ethers have positive liminf, and we believe that there are infinitely many.) 

As remarked earlier, many CA rules of interest provably exhibit more complex behavior for 
certain exceptional seeds. One important class of behavior is formalized by the following concept 
introduced in |GGP] . We call a seed ^o (or ^ configuration ^) a quasirepUcator with ether r] 
if the following holds. For some exceptional set of space-time points Q ^ A, each bounded 
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Figure 1.3: Four examples of Piggyback evolution: two replicators (with enlarged regions showing 
different ethers) from random seeds of length 30; a quasireplicator with seed 11111; and an 
apparently chaotic example with seed 100011011. 



1 INTRODUCTION 6 

component oi'I?\Q is filled with rj, while for some a > 1, the set a^^Q converges as n — )• oo in 
Hausdorff metric to a set of HausdorfF dimension strictly less than 2. 

Theorem 1.3 (Quasireplicators). For some diagonal- compliant and wide- compliant web CA 
rules, including Web-Xor and Piggyback, there exist seeds that are quasireplicators but not repli- 
cators. 

Examples of (provable) quasireplicators include the second and third seeds in Figure 11.21 
and the third seed in Figure 11.31 Certain other seeds appear to be neither replicators nor 
quasireplicators, but exhibit apparently chaotic behavior, although proving this seems very 
challenging. The fourth examples in each figure are in this category. In some very special cases 
we can prove chaotic behavior in a certain conditional sense, even for an infinite family of seeds 
whose number grows exponentially with their length. We discuss these issues further in the next 
section. 

Theorem 11.11 describes the space-time configuration away from A, and moreover states that 
this description is insensitive to 2s in the initial configuration. However, the result provides no 
information about the configuration close to A. The next result addresses this. The forward 
cone of a space-time point (x, t) is the set {(y, s) : \y — x\ < s — t}, and the forward cone of a 
set is the union of the forward cones of its points. 

Theorem 1.4 (Stability). Consider a diagonal- compliant or wide- compliant web CA, started 
from a uniformly random binary seed on [0,L]. With probability converging to 1 as L —^ cxo, the 
configuration of ^ in the forward cone of [0,-L] x {[ClogLj} is unchanged if any set of Os in ^o 
are changed to 2s. Here C is an absolute constant. If the CA has no spontaneous birth then 
with probability converging to 1 the same cone contains no 2s. 

We next discuss some ideas behind our proofs. Since in a web CA the web of Is evolves 
according to 1 Or 3, it easily follows that all Is lie in A(L). In the situation of Theorem 1 1.1 1 we will 
prove that immediately above each bounded component of Z^\A(L) there is a strip which blocks 
information flow. Furthermore, each such strip contains a spatially periodic configuration of Is, 
with the repeating unit being identical for all strips up to translation. This is a probabilistic 
statement, not a deterministic one, and the height of the strip is random. It will be proved 
using techniques of percolation theory. In contrast with classical percolation, the space-time 
configuration A of ^ Or 3 is not i.i.d., but has long-range dependence. We will make use of the 
key percolation result below, which we believe is interesting in its own right. 

A path is a finite or infinite sequence n of space-time points {xo,to), (xi, ti), . . . , {xn,tn) (,•••) 
with tj+i = tj + 1 and |xj+i — Xj| < 1 for all i. A path is diagonal if it satisfies |xj_|_i — Xi\ = 1 
for all i. Suppose Aq is given, and let A be the resulting configuration of 1 Or 3. We say that 
a path vr is empty if A(x, t) = for every (x, t) on tt. A path is ■wide if it is empty and it 
makes no diagonal step between two Is, i.e. it has no two consecutive points {x,t),{y,t + 1) 
with \x — y\ = 1 but X{y,t) = X{x,t + 1) = 1. (As suggested by the terminology, diagonal- 
compliance and wide-compliance of web CA refer to information fiow being restricted to paths 
of the appropriate type) . We now assume that the initial configuration Aq is uniformly random 
on Z, that is, Ao(x) takes values 0, 1 with equal probabilities independently for all x E Z. 
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Theorem 1.5 (Subcriticality). Consider the 1 Or 3 CA from a uniformly random initial con- 
figuration on 7L. We have 

(1.1) P(3 an empty diagonal path from Z x {0} to (0,t)) < e~^ , t > 0, 

for some absolute constant c > 0. The same conclusion holds for the existence of a wide path. 

In contrast, we prove that empty paths do percolate. 

Theorem 1.6 (Supercriticahty). For the 1 Or 3 CA from a uniformly random initial configu- 
ration on 7L, 

P(3 an infinite empty path from (0,0)) > 0. 

We now briefly discuss background to our resuhs. As remarked earher, CA that exhibit 
chaotic behavior for typical seeds but regular behavior for some seeds are apparently very com- 
mon. Empirical evidence strongly suggests that the one-dimensional rules Exactly 1 [GG3] . 
Perturbed Exactly 1 |GGPj . and EEED |GG4| are all in this category. It is natural to postulate 
a mechanism for this phenomenon, whereby chaos nucleates from certain local patterns, and, 
once started, invades all non-chaotic regions. It is tempting to conclude that this robustness of 
chaos might be universal law, akin to the second law of thermodynamics. 

To our knowledge the first compelling evidence to the contrary was presented in |GG2j . where 
a CA later called Extended 1 Or 3 was introduced. This rule arises naturally as the "2-layer 
extremal boundary dynamics" of a classical two-dimensional CA rule. Box 13. Piggyback is 
also the 2-layer extremal boundary dynamics of a two-dimensional rule. See Section [2] for more 
information. Extremal dynamics have been utilized very effectively in the analysis of Packard 
snowflake CA in |GGn[(TG2] . 

Extended 1 Or 3 was proved in |GGP] to admit both replicators and quasireplicators, and 
observed to generate apparent chaos from some seeds. Empirical evidence was presented that 
long random seeds are replicators with high probability, and thus that it is the ordered phase 
that is resistant to disorder. In this article we provide the first rigorous demonstration of this 
phenomenon. The classes of CA that we consider are strongly inspired by Extended 1 Or 3. 
We have not succeeded in proving that the conclusion of Theorem 11.11 holds in the case of 
Extended 1 Or 3, although this would follow if a certain natural conjecture (Conjecture 15. 2p 
were established. 

We note that the disorder-resistance phenomenon under consideration is somewhat reminis- 
cent of insensitivity of CA rules to random noise in the update rule, as in [GaclllGraj . 

Much CA research has focussed on evolution from carefully chosen initial configurations — 
a notable rigorous example is |Coo| . In contrast, rigorous results for CA from a random initial 
configurations are scarce, despite their potential importance in understanding self-organization. 
Most such research has been focused on nucleation, that is, random formation of centers that 
orchestrate a takeover of the available space. Notable examples include bootstrap percolation 
[Holl IBBDMj and excitable media models [FGG| . We also mention two previous works on 
additive dynamics started from a product measure, [Linj and |ENj : the latter finds an embedded 
random walk by an argument somewhat related to the methods in Section [5l 
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In many cases, percolation with long range dependence is extremely challenging to analyze 
rigorously (see [Wint iBBSj IGac2t IBSj . and references therein). Nevertheless, in our setting it 
turns out that the additivity of 1 Or 3 allows certain judiciously chosen percolation arguments 
to be carried through with relative ease. Translating results from an infinite random initial 
configuration to finite seeds also appears daunting, since the number of random bits is now finite. 
However, additivity introduces extensive periodicity and repetition into the configuration. With 
care, these properties can be used to advantage. This extreme form of long-range dependence 
provides the link between lack of percolation and evolution from random seeds, and is also the 
reason for formation of ethers. 

While our results provide a reasonably comprehensive picture of subcritical percolation be- 
havior for certain path types (diagonal and wide), it should be emphasized that the behavior 
for paths of supercitical type (empty paths) in the evolution from finite random seeds is not well 
understood. We discuss open questions and prove some preliminary results in this direction in 
Section [6l 

The article is organized as follows. In Section [2] we establish terminology, including the 
formal definitions of diagonal-compliance and wide-compliance, we introduce and discuss some 
further examples of CA having these properties, and we discuss how Theorem 11.31 is proved. 
Sections [3H3 are concerned entirely with properties of the additive rule 1 Or 3, from which 
properties of web CA are deduced later. In Section [3] we review properties (most of them well 
known) of i Or 3 started from a single occupied site, and in Section [H we use additivity to 
deduce basic properties of the evolution from random configurations. In Sections [5] and [6] we 
prove the percolation results. Theorems 11.51 and 11.61 respectively, and discuss other facts and 
open problems concerning percolation. In Section [7] we deduce key results about evolution of 

1 Or 3 from random seeds. Finally we return to web CA. In Section [8] we deduce Theorems ll.il 
and 11.41 and in Section [9] we prove Theorem 11.21 and show how to compute lower bounds on 
ether probabilities. 

2 Definitions, examples, and preliminary results 

2.1 Basic conventions 

Throughout the paper, A denotes the 1 Or 3 CA, while S, denotes a web CA. All our intervals 
will be subsets of Z or of Z x {t} for some t > 0. We adopt the convention that [a, 6] = and 
[a, b] X {t} = whenever b < a. 

Throughout, a site or a cell will refer to an element of Z; a point will be an element of 
space-time Z x [0,oo) C Z^. The state of a CA ^ at cell x and time t is denoted £,t{x) or 
^(x,t), depending on whether our focus is on time evolution or the space-time configuration. 
When specifying a seed, we always assume that all states left unspecified are 0. In diagrams of 
space-time evolution, state is colored white, state 1 is black or grey, and state 2 is a different 
non-greyscale color for each CA rule. 

We say that a collection of {0, l}-valued random variables is uniformly random if they are 
independent and take values and 1 each with probability 1/2. 
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2.2 Compliance 

In this section we formally introduce various families of web CA. As mentioned already, these 
will have 3 states and range 2. Thus the state of a site is ^t{x) S {0,1,2} for x G Z and 
t S [0, oo), and the evolution is given by 



^t+i{x) = f(Ct{x - 2),et(x - l),6(x),6(x + l),Ct{x + 2) 

for some function /. 

We reiterate our standing assumption that the Is of ^ behave as the 1 Or 3 CA. More 
precisely, writing 

(5(a) := l[a = 1] = amod2, a = 0,1,2, 

we assume that 

(2.1) <5(/(a, b, c, d, e)) = 5(6) + 6{c) + 6{d) mod 2 
for all a, b, c, d, e. Thus, if we define 

(2.2) Xt{x) := 5(6(x)), 

then (j2.ip implies that A satisfies the 1 Or 3 CA rule. We sometimes call A the first level of 
the process. We call a CA rule that satisfies (j2.ip and /(0, 0,0, 0,0) = a web rule. 

We now consider various further conditions that may be imposed on /. The idea will be 
that the flow of information concerning the distinction between states and 2 is blocked by 
Is (in various locations). Throughout the following, we take a,b,c,d,e and a' ,b',c',d' ,e' to be 
arbitrary satisfying 5{a) = 5{a'), 5{b) = 5{b'), etc. 

We say that the rule / is empty-compliant if 

/(a, 6, c, d, e) = f{a', b, c, d, e); 

that is, a cell's next state S,t+iix) depends on non-adjacent cells ^t(x±2) only through their first 
level. (Recall that by (j2.ip . the first level of the next state cannot depend on the non-adjacent 
cells at all). Similarly, we say that the rule is diagonal-compliant if 

/(a, b, c, d, e) = f{a', b, c, d, e'). 

It will be convenient to express the next conditions in terms of the new first-level states of 
the neighboring cells. Thus we denote 

£:=6{a)+6{b) + 6{c)mod2- 
r ■= 6{c) + 5{d) + 6{e) mod 2, 

so that if {a,b,c,d,e) = {^t{x-2),. . . ,^t{x + 2)) then (£,r) = {Xt+i{x -1), Xt+i{x + l)). We say 
that / is wide-compliant if it is empty-compliant and 

c = r = l implies f{a,b,c,d,e)=f{a',b,c,d',e'), 
and c = ^ = 1 implies f{a,b,c,d,e)=f{a,b',c,d,e). 



2 DEFINITIONS, EXAMPLES, AND PRELIMINARY RESULTS 10 



In a configuration A of i Or 3, a patli is said to be 9-free if it is empty and it contains no 
point (x, t) whose 5-point neighborhood {{x ± 1, t), (x it 1, t — 1), (x, t — 1)} contains 9 or more 
Is. Finally, we say a CA rule / is ^-free-compliant if it is empty-compliant and 

5{h) + 5{c) + 5{d) + 5{i) + 5(r) > 9 
implies /(a, b, c, d, e) = /(a', 6', c', d' , e). 

Recall the definition of no spontaneous birth from the introduction; this is equivalent to the 
condition that /(a, 6, c, d, e) / 2 whenever a, 6, c, d, e G {0, 1}. 

As suggested by the terminology, the behavior of cellular automata satisfying the above 
conditions is constrained by paths of the appropriate types. 

Lemma 2.1 (Compliance). Consider a web CA that is empty- compliant (respectively: diagonal- 
compliant, wide-compliant, or 9 -free- compliant). Consider two initial configurations £,o,£,q whose 
first levels agree (i.e. 6{^o{x)) = S{^q{x)) for all x), and define the first-level dynamics A via 
(j2.2p . Fix a point {y,t). If A has no empty path (respectively: empty diagonal, wide, or 9-free 
path) from any (x,0) at which ^0(3^) 7^ Co(^) ^^ {y^i)} then £,t{y) = Ctiv)- Moreover, if the CA 
has no spontaneous birth, then ^t{y) 7^ 2. 

Proof. Suppose, to the contrary, that £,t{y) 7^ Ct{y')- We need to show that there exists a path 
of the appropriate type from Z x {0} to (y, t). By induction, it suffices to exhibit the final step 
on this path. This is a straightforward verification. 

To prove the final claim in the no spontaneous birth case, consider the initial state .^g in 
which every 2 of .^0 is changed to 0. Then S,t{y) = C<(y) = 0- ^ 

Lemma 2.2 (3-free paths). In any configuration X of 1 Or 3, any 3-free path is wide. Any 
3- free- compliant web CA rule is wide-compliant. 

Proof. Assume that a 3-free path makes a leftward diagonal move on two space-time points in 
state 0. Denote the states a, b, c at nearby points thus: 

abO 
Oc 

We need to show that b and c cannot be both 1. However, if 5 = 1, then also a = 1, but then 
c = as the path is 3-free. This establishes the first claim. A similar argument gives the second 
claim. n 

We now state a simple but important lemma that says that, although the web rules have 
range 2, empty-compliance ensures that the "light speed" is essentially 1. 

Lemma 2.3 (Light speed). Assume an empty- compliant web CA. The state S,{x,t) depends on 
the initial configuration ^0 only through the states 

Xo{x-t-l), Co{x-t),...,Co{x + t), Xo{x + t + l), 

where X is defined by (|2.2p . 
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Proof. The given states determine the following states at time 1: 

Ai(x-t), Ci{x-t + l),...,^iix + t-l), Xi{x + t). 
Then we use induction. D 

2.3 Examples of rules 

We will introduce several examples of web CA, chosen to represent various behaviors. Finding 
such rules is not particularly difficult, and we know of many others with similar characteristics. 
Let the function A^i (resp. A'^2) count the number of Is (resp. 2s) among its arguments, and 

A^i2 = Ni + N2. 

Our first example is Web-Xor, whose update rule is given by 

'1 (b + c + d) mod 2 = 1 
f{a,b,c,d,e) = {2 (6 + c + d) mod2 = and N2ib,d) = 1 
otherwise. 

It is easy to check that Web-Xor is diagonal-compliant and has no spontaneous birth. Examples 
of its evolution are given in Figure [L2l The top example represents typical behavior: replication 
with zero ether from a long random seed. The middle two examples are quasireplicators, one very 
simple and one similar to the one in Theorem 8 of |GGP| . For many seeds including these two, 
quasireplication can be rigorously proved via inductive schemes that completely characterize the 
configuration at certain specified times. In more complicated cases, such schemes can be very 
laborious to construct, while in other cases it may be difficult even to determine whether the 
seed is a quasireplicator. We will not give proofs of quasireplication; instead we refer the reader 
to [GGP| for two typical examples of inductive schemes that feature in such arguments. We 
believe that the final example in Figure [L2] is chaotic. 

Even this simplest of rules displays a remarkable variety of behavior from "exceptional" 
seeds. Other interesting seeds that we have found include: 110010012 (a replicator with non- 
trivial pattern of 2s in the web), 110011112 (a quasireplicator with scale factor a = 4), 111001112 
(perhaps chaotic or a very complicated quasireplicator), 10110112 (apparent chaos restricted to 
one side). 

Modified Web-Xor also has no spontaneous birth, but the 2s obey a symmetric two-point Or 
rule in the presence of Is: 



f{a,b,c,d,e) 



1 {b + c + d) mod 2 = 1 

2 {b + c+d) mod 2 = 0, and 

either N2{b,d) = 1 or [N2{b,d) > 1 and Ni{£,b,c,d,r) > 1] 
otherwise. 



As seen in Figure [2711 this rule is capable of "mixed replication" with two different ethers (top). 
Note that Theorem 11.11 implies that with high probability this does not happen for long random 
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Figure 2.1: Modified Web-Xor with seeds 11111112 and 210001. 

seeds. The bottom example is apparently a quasireplicator, although we have no proof, and 
it seems that the inductive methods of |GGP| do not apply. Here and in the last example of 
Figure [L2| it is plausible that the evolution is driven by the advance of a front that lags behind 
the edge of the light cone by a power law. We will discuss this phenomenon in Section [6l 

In Web-adapted Rule 30, 2s evolve according to Rule 30 |Wol2] . except that 2s perform 
the three-point Or rule in the presence of Is when a neighborhood occupation number is small 
enough: 



f{a,b,c,d,e) 



1 {b + c + d) mod 2 = 1 

2 (6 + c + (i)mod2 = andiVi2(£,fe,c,d,r) <2, and 

either W3o[52{b),52{c),52id)] = 1 
or [N2{b,c,d) > 1 and Ni{£,b,c,d,r) > 1] 
otherwise. 



Here, w^o is the update rule for Rule 30, given by it;3o(ai, 02, ^3) = (^i + 02 + 03 + 0203) mod 2, 
and ^2(0) := l[o = 2]. Web-adapted Rule 30 is 3-free-compliant (and therefore wide-compliant) 
and has no spontaneous birth. See Figure [2^ for an example. One can prove that this instance 
is not a replicator, but is it chaotic? There are no known methods to prove chaotic evolution, or 
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Figure 2.2: Chaotic behavior of Web-adapted Rule 30 with seed 100010201. 



even universally agreed definitions of the concept; however, suppose one accepts the reasonable 
premise that Rule 30 generates a chaotic configuration p started from a single 1 [Wol2j . Then 
the example in Figure 12.21 is equally chaotic, in the sense that its evolution provably features 
larger and larger regions of p, at specific locations that are easily characterized. We will also 
show in Section [5] that an exponentially growing family of seeds exhibit conditional chaos in the 



same sense. 



The above rule may be modified in various ways so as to include spontaneous birth, resulting 
in further rules where Theorem 11.11 applies, yet in which many provable replicators have ethers 
with very long temporal period, perhaps too long to be seen experimentally. In the interests of 
brevity we omit the details. We briefly discuss bounds on the period in Section [8l 

A number of web rules arise naturally in analysis of two-dimensional CA, as we now explain. 
Consider a binary CA C,t £ {0, 1}^ , in which the new state of cell z is given by a rule defined 
on the Moore neighborhood M{z) := {z' ^ 1? : \\z' — z\\ao < 1}. We assume that state is 
quiescent, and that the CA solidifies, that is, Ct{z) = 1 implies C,t+i{z) = 1; the CA rule then 
only needs to specify when a z G Z^ becomes occupied, i.e., changes its state from at time t to 
1 at time t + 1. To each such CA we associate extremal boundary dynamics (EBB): assume 
that Co vanishes on Z x [1, oo) and let Aj be given by (^^ on Z x {t}. Observe that Aj is a one- 
dimensional CA whose space-time configuration is a lower bound on the final configuration 
Ccx) = Ut>oCi- Now assume that we extend the boundary layer to width 2, which leads to the CA 
it G {0, r,2}^ with the following rule: it{x) = 1 if C,t{x,t) = 1 (so that At = ^tmod2), it{x) = 2 
if Ct{x,t) = but (t+i{x,t) = 1, and Ct{x) = otherwise. Again, ^t is a one dimensional CA. As 
C,t{x,t — 1) = 1 exactly when either ^t-i{x) = 1 or it{x) = 2, ^j indeed determines two extremal 
layers of Ct and is thus called two-level EBD. The evolution of ^j also provides a lower bound 
on Qoo and is often useful when the bound provided by At "leaks" |GG2] . To conform with the 
rest of the paper, we assume throughout that the EBD is the 1 Or 3 CA. 

The natural setting for study of the issues addressed in this paper are general web CA, a 
much larger class than the two-level EBD rules. The latter, however, provide many interesting 
examples. In fact, the different ethers, quasireplicators and (apparent) chaotic behavior were 
first observed in the two- level EBD generated by the Box 13 solidification CA |GG2] . in which 



2 DEFINITIONS, EXAMPLES, AND PRELIMINARY RESULTS 14 

z becomes occupied at time t + 1 when the the number of occupied cells in M{z) at time t is 1 
or 3. The corresponding two-level EBD is called the Extended 1 Or 3 CA, and is given by 



f{a,b,c,d,e) 



1 {b + c + d) mod 2 = 1 

2 (6 + c + d)mod2 = and iVi2(^,r, 6,c, d) G {1,3} 
otherwise, 



as is easy to check; therefore, this rule is equivalent to the one with the same name introduced 
in |GGP] . This rule is 4- free-compliant, and is not covered by our main theorems. However, we 
establish some rigorous results in Sectional 

For simplicity, assume that the two-dimensional CA ^ is isotropic, that is, that its rule 
respects all isometries of the lattice Z^. Then there is a convenient sufficient condition that 
assures wide-compliance for its two-level EBD: when the neighborhood configuration is 

ale 
b 1 


the next state at the center cell is independent of c (i.e., depends only on a and b). This holds, 
for example, for the following solidification rule, which we call Perturbed Box 13. Given Q, let 
occi(2;) (resp. occoo(-z)) count the number of occupied cells among the four nearest neighbors 
of z (resp. in J\f{z)); then z becomes occupied if either 

• occi(z) = 2, or 

• occi(z) < 1 and occoo(-z) € {1,3}. 
See Figure 12.31 for an example. 

The resulting two-level EBD has the update rule 

' 1 {b + c + d) mod 2 = 1 

{b + c + d) mod 2 = 0, and 
/(a, 6, c, d,e) = { either Ni2{i, c, r) = 2 

or [Ni2{i,c,r) < 1 and Ni2{e,b,c,d,r) G {1,3}] 
otherwise. 

We call this web CA Piggyback. It is easy to see that it is wide-compliant, and has spontaneous 
birth. The top two examples in Figure [TTHl start from long random seeds and are replicators with 
diff'erent ethers. (We will have more to say about ethers for Piggyback in Section [9l) The third 
example is provably non-replicating, as it is a quasireplicator. The bottom example appears to 
be chaotic. Like the bottom picture in Figure 11.21 the evolution displays a tantalizing mixture 
of order and disorder. 

Our results on Piggyback have rigorous implications for the two-dimensional Perturbed Box 
13 rule (and similarly in other cases where 2-level EBD satisfies the conditions of Theorem II. ip . 
Here we summarize some initial observations, noting that further investigation is warranted. As 
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Figure 2.3: Perturbed Box 13 started from a seed in the square [0, 16]^. Initially occupied cells 
are black, and subsequently occupied cells are red or blue if they have state 1 or 2 respectively 
in the 2-level EBD, and otherwise grey. Unoccupied cells are white. 



suggested by Figure [231 the evolution of Perturbed Box 13 from a seed in [0,-L]^ is governed 
by four space-time configurations of Piggyback in four quadrants with boundaries at 45° to the 
axes. Depending on the behavior of each, we may make deductions about the final configuration 
Qao- In the case of a replicator with the "solid" ether (2)°°, as in the bottom quadrant in this 
example, clearly no further filling of the ether is possible after the second level of the EBD. By 
Theorem 11.21 it follows that that Perturbed Box 13 started from a uniformly random seed in 
[0, L]^ results in a final configuration Qo^ of density 1 in 1? with probability bounded away from 
as L — >■ oo. Certain other ethers of Piggyback can also be shown to fill in in a predictable 
manner, resulting in a corresponding ether for Perturbed Box 13, as in the top quadrant. A 
similar analysis can likely be carried through for certain simple quasireplicators such as the one 
in the right quadrant. When Piggyback is a replicator with zero ether, as in the left quadrant, 
it appears plausible that the subsequent filling- in by Perturbed Box 13 results in a chaotic final 
configuration. See |GGll[GG2j for detailed analysis of the filling-in process for some other EBD. 

We conclude by mentioning a natural rule that seems intractable by our current methods. 
Web 1 Or 3 is the web CA in which 2s perform 1 Or 3 on the points not occupied by Is: 



f{a,b,c,d,e) 



'\ (b + c + d) mod 2 
2 (b + c + d) mod 2 
otherwise. 



1 

and N2{b, c, d) mod 2 = 1 
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Figure 2.4: Chaotic behavior of Web 1 Or 3 from a random seed of 32 sites. 

Figure [23] gives an example of an evolution from a random seed of Is and 2s, with an apparent 
message of near-criticahty and chaos. 

2.4 Generalizations 

The simplest additive rule, Xor CA jit, is defined on the state space {0, 1}^ by 

fitix) = iJLt-i{x) + P^t+i{x) mod 2. 

One might consider /i, and not A, to be the most natural candidate for the web dynamics. 
However, while // does have some points of interest (see, for example Proposition 16. 5p . many of 
the main issues we consider become trivial in this setting. For example, // either only occupies 
points satisfying a parity constraint or generates an impenetrable web even for empty paths 
[BDRllGGl] . 

In the other direction one might ask whether similar results hold if A is replaced by an 
arbitrary additive rule. It is indeed likely that a more general theory could be developed in this 
setting. One complication is that predecessors of the all-0 state will no longer necessarily be 
unique (as they are for 1 Or 3 — see Lemma [3. 4p and as a result "mixed replicators" similar to 
the top example of Figure 12.11 may be the norm. 

On the other hand, all our results generalize with appropriate minor changes in the definitions 
to CA with a quiescent state 0, first-level state 1 and other states 2, . . . , s. 

3 Additive dynamics from a single occupied site 

Recall that A* denotes 1 Or 3 started from a single 1. In this section we collect properties that 
we will need. All these results are elementary and many are well known. First is a rescaling 
property, illustrated in Figure ISTTl 
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Figure 3.1: An illustration of Lemma [3 -H with m = 2. Highlighted points comprise a "separated 
out" copy of A*. 

Lemma 3.1 (Rescaling). For any nonnegative integers a and m, 

,. _\K{y) ifx = 2"^yforyeZ 
I otherwise. 

Proof. The case m = 1 follows from additivity on observing that A* is 10101. For m > 1 we 
apply the m = 1 case iteratively. D 

Lemma 3.2 (Periodicity properties). 

(i) For t > 0, Ai*(0) = Xt{±t) = 1 while Xti±{t - 1)) = tmod2. 

(ii) For all n >0, A2n(x) = 1 exactly at x = 0, ±2". 
(iii) For all n>0, Xl^_^_^„_i{x) = 1 exactly at x = 0,±2",±(2" + 2"-^). 

(iv) For any k > 1, the sequence of edge configurations of X* on [t — k + l,t] x {t} is periodic 
(from f = on) with period equal to IP where IP"^ <k <'2P . 

Proof. Parts (ii) and (iii) follow from Lemma [3 -H and (iv) follows from (ii), with (i) as a special 
case. D 

For some purposes, the following recursive description of A* is useful, a variant of the one 
given |Wil] . See Figure [331 for an illustration. Given a space-time configuration ^ on 5„ = 
[0, 2"] X [0, 2" — 1], we say that A is placed at a space-time point s if the configuration in s + Sn 
is the corresponding translate of A. Let Bn be the space-time configuration of A* on Sn. Reflect 
Bn around its vertical bisector and denote the resulting configuration on Sn by -B„. 

Lemma 3.3 (Recursion). We have Bq = 10 and Bi = -.-.„• Moreover, for n > 2, Bn is 

obtained by placing B^-i at (0,0) and at (2"-i,2"-i); B„_2 at (0,2"-^) and at (0, 2""! -F2"-2); 
and Bn-2 at (2"'~^, 2*^"^) and at (2"~^,2"~^ -|- 2^~^). All placements result in consistent state 
assignments at overlaps. 
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Figure 3.2: Evolution of A* with highlighted boundary strip of width 8 and temporal period 8. 



I 



■ ■ ■!■ 

! ■ "B" ■ !b 



Figure 3.3: Recursive description of A*: B^ is composed of two copies of B^ (red), 2 copies of 
B2 (green) and two copies of B2 (blue). 



Proof. This follows easily from (i) , (ii) and (iii) of Lemma 13.21 



D 



Our results for seeds depend on the fact that A* has certain a unique periodic configuration 
above every region of Os. This property does not hold for general additive rules. 

Lemma 3.4 (Predecessors of 0). For an arbitrary initial state Xq, suppose that Xt = on [a, b], 
but Xt-i ^ on [a— 1, 5 + 1]. Then Xt-i is a subword of the periodic word (110)°° on [a — 1, 6+1]. 

Proof. Consider the 4 possible values for the pair Xt-i{a — 1) and Xt-i{a). Once these states 
are fixed, the rest of At_i on [a — 1, 6 + 1] can be determined sequentially. D 



Fix an initial state for A. A void is a finite inverted triangle of the form Uj>o([a + i, 6 — i] x 
{t + i}), on which the configuration is identically 0, and that is maximal with these properties 
with respect to inclusion. Its width is 6 — a + 1, and its start time is t. 
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Figure 3.4: Illustration of Proposition 13.61 with m = 2. The highlighted intervals at distance 2^ 
above two selected voids have the claimed periodic configuration. 

Lemma 3.5 (Voids). In A*, each void has width 2—1 and start time divisible by 2 for some 
integer k > 1. Furthermore, for every fixed k, the union of all voids of width at least 2^ — 1 has 
density 1 within the forward cone o/(0,0). 

Proof. This is a straightforward application of Lemma |3.3[ D 

Finally, we deduce the following fact, which will be crucial for our results on percolation and 
ethers. See Figure [33] for an illustration. 

Proposition 3.6 (Periodic interval above a void). In A*, assume that [a, 6] x {t} is the top row 
of a void of width 2—1. For m < k, the state of interval [a — 2"^, 6 + 2™] x {t — 2"^} is a segment 
of the following infinite periodic string of period 3-2™'.- 

(3.1) (1 dl 1 dl IIZI)°°. 

Here, I I represents a string of 2"^ — 1 consecutive Os, and the segment begins and ends with a 
full CH. 

Proof. As t is divisible by 2 , and therefore by 2*", A* on [a, b] x {t} is of the form 

dl CH ... dl, 

by Lemma 13.11 Then, by the same lemma, and Lemma 13.41 applied to Xt/2'^ ; the configuration 
on [a — 2™, 6 + 2™] x {t — 2*"} is either of the claimed type started and ended with CZI, or all 
Os. The latter possibility contradicts maximality of the original void. D 
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4 Duality and randomness 

When the initial configuration of i Or 3 is uniformly random (on some set), the resulting space- 
time configuration is of course not uniformly random but has a high degree of dependence. 
Nevertheless, in this section we show how to identify space-time sets on which the randomness 
is uniform. The additive structure of the CA rule ensures that the space-time configuration is a 
linear function (modulo 2) of the initial states, and the idea is to find cases where the associated 
matrix is upper triangular. 

Recall that A* is the 1 Or 3 rule started with only the origin occupied. Let A^ denote the 
rule started with the set of initially occupied sites exactly equal to A C Z. We will extensively 
use the following version of cancellative duality. 

Lemma 4.1 (Duality). We have Xf{x) = Yly^A ■^'(^ ~ v) ™od2. 

Proof. This follows easily by additivity and induction on t. D 

Observe that by symmetry and translation-invariance, A*(x — y) = A*(y — x) = Aj (y). 

Suppose we have an ordered set S = {{xi,ti) : i = 1,2, ...,n}, of space-time points. A 
function F : 5" — >■ Z is a dual assignment for S if for alH, j G {1, . . . , n}, 

X'{xj - F{xi,ti),tj) = 

(There is no restriction when j > i). We think of F{-, •) as sites in the initial configuration. The 
idea is that in order to determine A (xj, tj), we need new information about A at each successive 
i. 

Proposition 4.2 (Randomness via dual assignment). Suppose that the initial configuration Aq 
of 1 Or 3 is uniformly random on some fixed set K C.'L and deterministic on K^ . Let S be a 
fixed set of space-time points. If S has a dual assignment whose image is contained in K , then 
A is uniformly random on S. 




Proof. Writing 



and 



we have by Lemma l4.H 



Ki = {y e K : Xt^{x, - y) = 1}, 
Kl = {yeK'':XtM-y) = l}. 

Cj = ^ Ao(y)mod2, 
yeKi 



X{xi,ti) = ^ Ao(y) + Cimod2. 
y&Ki 

But Ki contains an element, F{xi,ti), that is not in IJj^iKj, therefore X{xi,ti) is uniformly 
random conditional on {X{xj,tj) : j <i). D 
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A particularly useful special case is that a 1 adjacent to a string of Os in A* heralds uniformly 
random intervals in the evolution from a random seed. 

Corollary 4.3 (Random intervals). Fix integers a and L,k > 0. Let the initial configuration 
Xq be a uniformly random binary seed on [0,L], and suppose that A* on [a,a + k] is 1 followed by 
k Os. Then for any x G [—a, L — a — k], the configuration Xt is uniformly random on [x, x + k]. 

Proof. To find a dual assignment of [x,x + k] x {t}, order the set from left to right, and let 
F{y,t) = y + a. Clearly, the image of this assignment is contained in [0,L]. Now apply Propo- 
sition 1121 □ 



5 Subcritical percolation 

In this section we prove Theorem ll.5l which states that when 1 Or 3 is started from a uniformly 
random initial configuration on Z, the probability of an empty diagonal or wide path from the 
initial interval Z x {0} to the point (0,t) decays exponentially in t. See Figures [5TT] and [5T3l for 
the diagonal and wide cases respectively. Note the contrast with Figure 16.11 in the next section 
for empty paths. 

Our approach is to use dual assignments to control the probabilities of paths, but the details 
of the argument are very different for the two types of path. A diagonal path has 2 choices at 
each step, and any given point has state with probability 1/2, suggesting a critical bound. To 
improve this to a subcritical bound we consider a leftmost path, and use special properties of 
A. On the other hand, we control wide paths via a random process of space-time intervals that 
terminates when an interval has even length. 

Later in the section we also discuss 9-iiee paths, and show that notwithstanding Theorem 1 1.5 1 
there is an exponential family of initial configurations for which percolation by wide paths does 
occur. 

5.1 Empty diagonal paths 

Proof of Theorem \1.5\ case of empty diagonal paths. We may assume without loss of generality 
that t is even, since the configuration at time 1 is also uniformly random, and the probability 
in question is strictly less than 1 for t = 1. 

Fix a diagonal path it from (x, 0) to (0,t). We will find an upper bound for the probability 
that TT is the leftmost empty diagonal path from Z x {0} to (0, t). To this end, partition the steps 
of vr into segments of length 2. During each such segment, the path has one of the following 
forms: left-left, right-right, left-right, or right-left. When it makes a right-left move, that is 
(x, s) — > (x + 1, s + 1) — >■ (x, s + 2), the leftmost property requires a 1 at (x — 1, s + 1); we call 
these points (which are not on the path) the corner points of the path, and let N(tt) be their 
number, i.e., the number of right-left segments that start at even times. 
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Figure 5.1: All empty diagonal paths from an interval at time are highlighted in blue. The 
initial configuration is uniformly random. 



-k 1 1 1 1 1 Ct 1 i — I — I — I — J — I — 7I- 

^ \ \ \ / 



Figure 5.2: A leftmost diagonal path together with a dual assignment for the points of the path 
(white discs) and the corner points (black discs). The dashed lines connect each point to its 
assigned position in the initial configuration. 



We will give a dual assignment F of the path together with the set of its corner points (see 
Figure [5^ for an illustration). Order points on the path with increasing time, and place a corner 
point (x, s) in the ordering immediately after the point on the path at the time s + 1. For every 
corner point (3;,s), let F{x,s) = x — s + 1. For every point {x,s) on the path, let F{x,s) be 
either x — s or x + s, according to whether the path arrives to {x, s) from the right (i.e., from 
(x + 1, s — 1)) or from the left (i.e., from (x — 1, s — 1)), respectively. We let F{x, 0) = x. 

To check that F is a dual assignment, we will use Lemma I3.2l fi). Fix a point (x,s) on 
the path. All positions assigned by F to points earlier in the order lie in [x — s + 2,x + s] or 
[x — s,x + s — 2] according to whether the path arrives to {x, s) from the right of left, so the 
required condition is satisfied for this point. Now suppose (x — 1, s + 1) is a corner point arising 
from the moves (x, s) — )■ (x + 1, s + 1) — )• (x, s + 2) in the path. This corner point is assigned to 
-F(x — 1, s+1) = X — s — 1. We have F{x, s+2) = x — s — 2, and all points earlier than (x, s+2) were 
assigned integers at least x — s. Since s + 1 is odd, A*(x — 1 — (x — s — 1), s + 1) = A*(s, s + 1) = 1. 
Finally, since s + 2 is even, A*(x — (x — s — 1), s + 2) = X*(s + 1, s + 2) = 0, as required. 
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Figure 5.3: All wide paths from an interval at time are highlighted in blue. The initial 
configuration is uniformly random. 



Now, using Proposition 14. 2| 

P(an empty diagonal path from Z x {0} to (0, t) exists) 
< y P(7r is the leftmost empty diagonal path from Z x {0} to (0,t)) 

where both sums are over all diagonal paths from Z x {0} to (0, t). Let Pt be the last sum above. 
Then, by considering the last two steps of the path, 

SO, recalling that t is even, Pt = (1/2) • (7/8)*/^. D 

As an aside, we mention that the assertion of Theorem 11.51 for diagonal paths also holds 
when A is replaced by the Xor CA //, with a much simpler proof, since the set of all space-time 
points that the origin is connected to by diagonal paths is a rectangle. 

5.2 Wide paths 

Proof of Theorem \1.5\ case of wide paths. We will prove that 

(5.1) P(3 a wide path from (0,0) to Z x {t} \ Ao(0) = O) < e'"* 

for some c > 0. This clearly suffices by translation-invariance, since there are only 2t + 1 points 
at time from which a path can reach (0, t). Therefore, we will henceforth assume that Ao(0) = 
and that Aq is uniformly random elsewhere. 

We recursively define intervals It = [Lt, Rt] for t = —1,0, ... ,T, where T < oo, as follows. 
Start with L_i = i?_i = 0; then let Iq be the maximal subinterval of Z containing on which 
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Figure 5.4: The process of intervals of zeros used to prove non-percolation by wide paths. Here 
the interval lengths (|Io|, |/i|, . . . , |/t|) are (5,3,1, 7,5,3,1, 4,2,0). The witness points are 
highlighted. The corresponding binary sequence (Xi, . . . ,X2t+i) is shown below; it is obtained 
by reading the states of the witness points in conventional text order on the page, except with 
the left (red) intervals reversed. The dual assignment of witness points to initial positions is 
shown via dashed lines (witness points in the top row are assigned to themselves). 

Ao = 0. If It = 0, then we set T = t and there is no It+i- Otherwise, if \It\ > 2, then It+i is 
the interval [Lt + l,Rt — 1] (which is when \It\ = 2). If |/t| = 1, then It+i is the maximal 
subinterval of Z containing Rt = Lt on which Xt+i = 0. Observe that for each t < T we have 
At = on /(, while Xt{Lt — 1) = \t{Rt + 1) = 1- This follows from the CA rule for \t by induction 
on t; the key observation is that if \It\ = 1 then Ai+i(Lt) = (see Figure [57i|) . Furthermore, any 
wide path started at (0,0) is within Ut^rih x {*})• 

We now define an ordered sequence of 2r-|-l space-time points, which we call ■witness points, 
associated with the above sequence of intervals. If \It\ = 1, we call i -|- 1 a refresh time; we 
also declare a refresh time. Let = tq < ti < . . . be the refresh times. We build the sequence 
of witness points by appending certain points at each refresh time r,, in order. Specifically, for 
every i, we append all points in I^-- x {tj}, with the exception of zi = (L^-.-i, tj), in the following 
order: points to the left of Zi in the right-to-left order, followed by points to the right of Zi in 
the left-to-right order. Let Xi = A(sj), where si, . . . , S2t+i are the witness points in the order 
described. Write X for the random finite or infinite sequence given by X = (Xi, . . . ,X2t+i) if 
T < oo and X = (Xi, X2, . . .) if T = 00. Our goal is to show that X is equal in distribution to 
a sequence of independent fair coin flips stopped at a certain a.s. finite stopping time. 

Let yi, 12, . . . be independent random variables taking values and 1 with equal probability. 
Partition this sequence into blocks of the form 0'*10''1, where a,h > 0, and let S" > 1 be the 
location of the endpoint of the first such block of odd length. Then S is a.s. finite and 



(5.2) P(5>t)<e 



-ct 
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Figure 5.5: The sequence of principal voids 14, V2, . . . , Vg, numbered from top to bottom; Vi 
and V2 each consist of a single point. 

since S is at most the waiting time for the pattern 11011. Write Y' := (Yi, . . . , l^). We claim 
(5.3) X = Y'. 

Once (j5.3p is proved, the exponential bound (|5.2p implies (j5.ip . since 2T + 1 = 5. 

We now proceed to prove (j5.3p . Since P(S' < cxo) = 1, the random variable Y' has countable 
support, thus it suffices to show that P(X = y) = P{Y' = y) for any y with P(Y' = y) > 0. 
Choose such a y. The event {X = y} determines T and Iq, ■ ■ ■ ,It, and therefore the locations 
of the witness points. It follows that the event {X = y} is precisely the event that A takes the 
specified values y on these (deterministic) witness points. Now we use a dual assignment to 
show via Proposition 14.21 

/ 1 \ length of y 

PiX = y)= (-) 

as required. Considering the witness points in their order, we assign to each (x, t) either x + t 
or X — t, depending on whether it is to the right or left of Zi in its interval. (See Figure [531) 
This is a dual assignment simply because A*(±i) = 1 and A*(a;) = for x G [—t,t]'^. D 

5.3 ^-free paths 

Since 3- free paths are wide (Lemma 12. 2p . the exponential bound (jl.ip holds for 3-free paths, 
and there are no infinite 3-free paths when Aq is uniformly random on Z. Do such paths exist 
started from special initial conditions? Indeed they do, as shown by our next result. 

Define the sequence of principal voids Vi,V2, ■ ■ ■ of A* according to Figure 15.51 and for 

L > let Wi = Viri{Vi + (L, 0)) (the notation means that Vi is translated by the vector (L, 0)). 

Proposition 5.1 (Exceptional percolation). Assume Aq is a seed on [0,-L] which vanishes out- 
side 4Z, where L is a multiple of ^. Define the sets Wi as above. Let i be such that 2* > 16 L^ . 
From every point in Wi there is a 3-free path to some point in Wj+i. In particular, from any 
such point there is an infinite 3-free path. 
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Proof. For a maximal interval If C Z on which A^ vanishes, define its successor It-^-i to be the 
maximal interval on which A^+i vanishes and It n It+i ^ 0; if such an interval does not exit, let 
It+i = 0. 

Observe that At vanishes outside 2Z at all even times t. Using this, it is easy to verify by 
case-checking that every maximal interval of Os has odd length at every time. It follows that any 
nonempty maximal interval of Os has a nonempty successor. Similarly, since Xt vanishes outside 
4Z when t a multiple of 4, it is easily verified that for any t, from every point in It x {f } there 
is a 3-free path to some point in It+i x {t + 1}. 

Let the apex Wi be the (unique) bottommost point of Wi. Clearly, from every point in Wi 
there is a 3-free path to Wi. Moreover, by the above, there is a 3-free path from lOj to Z x {ipi/i+i }> 
where twi+i is the time of the top interval of VFj+i. The condition 2* > 16L^ ensures that this 
top interval contains the intersection of the forward cone of Wi with Z x {tvFi+i}- D 

The above result implies furthermore that there are, for some c > 1, at least c seeds on 
[0, L] for which there is a 3-free path from Z x {0} going through all Wi for i sufficiently large. 
This follows by replacing a seed Aq satisfying the condition of the theorem with its successor Aj 
for t a sufficiently large multiple of L. If all Os within such a seed are replaced by 2s, this also 
yields an exponential family of seeds for the Web-adapted Rule 30 CA with the same chaotic 
properties as discussed in Section [21 

We conclude this section with the following conjecture supported by computer experiments. 

Conjecture 5.2. The exponential hound (II. ID in Theorem I J.5I holds when diagonal path is 
replaced by ^-/ree path. 



6 Supercritical percolation 

In this section we consider empty paths from random initial conditions, and in particular we 
prove the percolation result Theorem 1 1.61 The results of this section are not needed for the proofs 
of Theorems ll.lfll.41 but they are of independent interest and complement those of the previous 
section. We also consider initial conditions where the randomness is restricted to the half-line 
or a finite seed. Here many questions are open, but we establish some preliminary results. The 
questions we consider are relevant to further understanding certain web CA behavior. 

6.1 Percolation of empty paths 

As Figure 16.11 suggests, the set of points reachable by empty paths emanating from an interval 
at time form an interval at each subsequent time. With random initial conditions, this interval 
spreads linearly provided it survives. Proving this is the key to Theorem 11.61 

Suppose Aq is given. The rightward Z-path from a space-time point (x, t) is an infinite 
sequence of points (^s, s), s > t defined as follows. Start with rt = x. Inductively, let r^+i be the 
largest integer y in (— cxo, r^ + 1] for which Xs+i{y) is 0; or if there is no such y we take r^ = — oo 
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Figure 6.1: All empty paths from an interval at time are highlighted in blue. 

for all u > s. Note that X(rs,s) = for all s > t for which r^ is finite, but not necessarily 
for s = t. Analogously, we define the left'ward Z-path {ls,s), s > t hy reversing the space 
coordinate in the definition. 

Lemma 6.1 (Properties of Z-paths). Suppose X is 1 Or 3 from any initial configuration. 

(i) Suppose A(0, 0) = and let {rt,t), t > be the rightward Z -path from (0,0). If x < r^ and 
X{x,t) = then there is a empty path from (— oo,0] x {0} to {x,t). 

(ii) Fix an interval [a,b] with a < b. Let (^i,t) be the leftward Z-path from, {a,0), and {rt,t) 
the rightward Z-path from, (6,0). Suppose that ig < r^ for every s < t. Then for any 
y £ [it, rt] with Xt{y) = 0, there is an empty path from [a — 2,b + 2] x {0} to (y, t). 

(iii) Under the assumptions o/(ii), suppose also that Ao(a) = Ao(6) = 0. Then for any y € [it, ft] 
with Xt{y) = 0, there is an empty path from [a,b] x {0} to {y,t). 

(iv) Conversely, if there is an empty path from [a,b] x {0} to some {y,t), then ig < r^ for all 
s < t, and it 1^ y 1^ ft- 

Proof. We omit the proof of (i) , as it is similar to the proof of (iii) , which proceeds by induction 
as follows. The argument reduces to verifying (iii) at time t = 1. Assume ii < ri and take 
y G [^ij^^i] C [a — 1,6 + 1] with Ai(y) = 0. Then there exists an a; G {y — l,y,y + 1} with 
Ao(x) = 0. It remains to verify that x can be chosen to be in [a, b]. If y G [a + 1, 6 — 1] this is 
clear; if y G {6, 6 + 1} we may take x = b and if y G {a, a + 1} we may take x = a. 

The above argument also proves (ii): we verify the claim at time t = 1 and then use (iii). 
The last claim (iv) is an easy consequence of definitions of empty and .Z^-paths. D 

The key fact in establishing percolation of empty paths is that rt has drift 1/4. The proof is 
somewhat similar to that of non-percolation for wide paths. Theorem 11.51 

Lemma 6.2 (Drift). Suppose that the initial configuration Aq is uniformly random on 'L\ {0} 
and Ao(0) = 0. Let {rt,t), t > be the rightward Z-path from (0, 0). For every e > there exists 
a constant c = c(e) > so that P{\rt — t/4\ > et) < e~^*. 
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Figure 6.2: The rightward Z-path (solid lines) from the origin, together with its refresh points 
(circled), and witness points (highlighted in red). A dual assignment of the witness points to 
initial positions is indicated by the dashed lines. The states of the witness points in the order 
they are examined are shown below. 



Proof. We first describe an exploration process that determines the rightward Z-path from the 
origin (0,0). We designate (0,0) to be the first refresh point. Now we examine the states 
of the points (1, 1), (0, 1), (— 1, 1), (— 2, 1), . . ., in this order, until we find the first point with 
state 0. Let G be the number of points examined, and call them witness points. Since the 
states of these witness points are 01 • • • 1 (from left to right), certain states at the immediately 
following time steps are determined. Specifically, the pattern 01 • • • 1 is immediately followed by 
patterns of the saine form, but with the length decreasing by 2 at each step and centered at 
the same location, ending with either 01 or according to whether G was even or odd. (See 
Figure [6^2]) . We designate the location of the in this last pattern to be the next refresh point. 
It is (1— [G/2J, \G/2]). Now iterate the process starting at the new refresh point. Note that the 
rightward Z-path from (0, 0) consists precisely of the Os at the left ends of the 01 • • • 1 patterns, 
including the refresh points. Observe also that the Z-path is determined by the locations of the 
refresh points, and that these are determined by examination of the witness points. 

Now consider the above exploration process for the initial configuration that is uniformly 
random on Z \ {0} and at 0. Let (Xj)j>i be the sequence of states of the witness points, 
in the order that they are examined by the exploration process. We claim that (Xj)j>i is 
uniformly random. It suffices to check that (Xi,...,X„) is uniformly random. This follows 
from Proposition 14. 2( by the dual assignment in which a witness point (x, t) is assigned to x + t 
if it is the rightmost witness point in its 01 • • • 1 pattern, and otherwise to x — t. See Figure [ 



Let Gi be the number of witness points examined in the row immediately immediately below 
the ith refresh point. Then (Gi) are i.i.d. Geometric(l/2) random variables. Furthermore, 
the sequence of refresh points is a random walk on Z^ with steps (1 — [Gj/2j, \Gi/2]). As 
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E[G,j/2j = 2/3 and E[Gj/2] = 4/3, each step has expectation vector (1/3,4/3). The proof is 
concluded by standard large deviation estimates. D 

Proof of Theorem MM For L to be chosen later, consider the leftward path (£i,t) started at 
(— L, 0) and the rightward path (r^, f) started at (0, L). Then, by a union bound and symmetry, 

P(^t < n Vt) > P(£i < -1 and n > 1 Vi) > 1 - 2P(ri < for some t). 

By Lemma 16.21 for L large enough we have P(r( < for some t) < 1/3. Call a site x € Z good 
if an infinite empty path starts at (x,0). Thus, by Lemma iG.lT iil. 

P([-L - 2, L + 2] contains some good site) >P{it< n Vt) > 1/3. 

Consequently, by translation-invariance, P(0 is good) > 1/[3(2L + 5)]. D 

6.2 Empty paths for half-lines and seeds 

How do empty paths behave when the initial configuration is a random seed? This question is 
largely unresolved. (In contrast, the next section will provide detailed answers for diagonal and 
wide paths). A first step would be to understand the case of a uniformly random half-line, for 
which the following conjecture is natural given Lemma l6.2i 

Conjecture 6.3. Suppose the initial condition Aq is uniformly random on [l,oo) and else- 
where. Let {rt, t), t >0 be the rightward Z-path from (0, 0). Then rt/t — > 1/4 as t ^ oo. 

We prove that a much weaker statement holds deterministically: for an initial configuration 
supported in a half-line, empty paths penetrate arbitrarily far into its forward cone. 

Lemma 6.4 (Unbounded penetration). Assume that the initial condition Aq has no Is outside 
[l,oo). Let {rt,t) be the rightward Z-path from (0,0). Then sup^(ri + t) = oo. 

Proof. We first observe that for any initial configuration Aq of 1 Or 3, if {x, t) has state and 
t > 1 then at least one of the three points {x,t — 1), (x it l,t — 1) has state also. Iterating 
this we see that there must be an empty path from Z x {0} to (x,t). We call any such path an 
ancestral path of {x,t). 

Now, under the conditions or the lemma, note that for any m >0, the sequence of configura- 
tions on the intervals [—t, —t-\-m + l] x {t} is periodic in t starting from some time tp depending 
on m and Aq- For a > 0, define the leftward diagonal Da ■= {{a — t,t) : t > 0}. Then A cannot 
be identically 1 on two consecutive diagonals Da and Da+i, and also cannot be identically 1 
on Da and identically on Da+i- (Indeed, in either case we deduce that A is also 1 on Da-i, 
leading to a contradiction by induction.) 

Fix m > 0. We will show that for some t there an empty path from (— oo,0] x {0} to 
{(— t + 771 + 1, t), (— t + m, t)}, which suffices by Lemma [6. 11 To verify this claim, we may assume 
that the periodic orbit commences initially, i.e., that tp = 0. There must be a time t with either 
Xt{—t + m + 1) = or At(— t + m) = 0; by periodicity there must be infinitely many such times. 
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Now take the leftmost ancestral path of one of these two points. Suppose this path does not 
start on (—00,0] x {0}. Then, if t is large enough, the path has a diagonal segment longer than 
the period of the orbit; additionally, all states immediately to the left of such a segment must 
be 1. By periodicity, we have, for some a G [0,m + 1], infinite diagonals Da and -Dq+i on which 
A is identically 1 and 0, respectively. This is in contradiction with our observations above. D 

We remark that the supremum in the above lemma cannot be replaced with a limit; a 
counterexample is Aq ^ 1 on [2,oo) and Ao(l) = 0. 

Returning to our earlier question on seeds, Figure 16.31 (top) shows the set of all points on 
empty paths from (—1, 0), when Aq is a random seed on [0, 25]. We believe that for typical long 
seeds, the right frontier of this set lags behind the right edge of the forward cone of the seed by 
a non-trivial power of t in the limit t — )■ oo. This is a natural guess, since the frontier has speed 
1 in the voids, but presumably speed 1/4 on the fractal set occupied by randomness. It appears 
plausible that such a process is a driving force behind the evolution of some exceptional seeds 
for web CA, including the examples in Figures 11.21 (bottom). [2T] (bottom), and possiblv 12.41 

We believe that similar power law behavior holds for some specific small seeds. One example 
is shown in Figure [6]3] (middle): the seed is 1000 0001, and empty paths from the middle 
are highlighted. However, some seeds exhibit entirely different behavior. The bottom picture 
shows the empty paths from the two Os in the seed 1000 00010000 0001. Despite apparent 
initial similarity to the previous case, here the rightmost point (rj,t) reachable at time t has 
rt/t bounded strictly between and 1 at t ^- oo. Indeed, the rescaled path 2~"'{(rt,t) : t > 0} 
converges as n — )• oo to a variant of the Cantor function or "devil's staircase". This may be 
proved by an inductive scheme. 

As a preliminary step towards the power law behavior postulated above, we prove a version 
in a simplified setting. Recall from the Section [2] that /i denotes the Xor additive cellular 
automaton rule. Given a configuration ^ G {0, i}^x|o.°o)^ ^e define the x-P^th starting from a 
point {x, 0) to be the sequence of points {{xt, t) : t >0) given by xq = x and 

jxt, ^i{xt,t) = l; 

Xf-^-l = < 

[xt + l, iJ.{xt,t) = 0. 

In other words, the path makes a down step from a 1, and a diagonal step from a 0. This is 
intended as a simplified model for a rightward Z-path, which moves with speed 1 in Os, but with 
a slower speed in a random configuration. 

Proposition 6.5 (Power law for Xor). Let /i be the Xor CA with initial configuration hq equal 
to 1 on the two-point set {—1,0} and elsewhere. The x-po-th {{xt,t))t>o starting from (0,0) 
satisfies 

Xt = t-e(ti°e2/iog3^ ast^oo. 

Proof. We first note some easy facts about /x. Denote the interval of points R{k,t) := {{i,t) : 
t — 2^^ < i < t) on the right side of the forward cone of the origin. For any A; > 1, the state- 
vectors {ii{z) : z € R{k,t)) on these intervals form a periodic sequence in t with period 2^^^^. 
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Figure 6.3: The set of all points (blue) on empty paths starting from certain initial points, in 
1 Or 3 started from three different seeds: two apparent power-law cases, and a devil's staircase. 
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Figure 6.4: The x-P^'th from the origin in the Xor ceUular automaton, together with the con- 
struction used in its analysis. The strips 5(2) (pink) and S{3) \ S{2) (orange) are shaded. 

Furthermore, the sequence of state- vectors on the intervals R{k + l,t)\ R{k, t) consists precisely 
of the all-0 vector repeated 2 times followed by the first 2 state- vectors for R{k,t) (all 
repeated with period 2^). See Figure [63] for an illustration of the case k = 2. 

Let -Efc := min{t > : {xt,t) ^ R{k,t)}; this is the time at which the x-path leaves the 
diagonal strip S{k) := [JtR{k,t). This can only happen at a down step, which can occur only 
at a 1 of /i in the leftmost diagonal of S{k). It follows that Ej. is divisible by 2^^^; write 
Ek = 2 Cfc. For example (referring to Figure [63]), we have 62 = 3 and 63 = 4. 

In order to leave the strip ^(A; + 1), the path must first leave S{k), and then leave S{k + 1) \ 
S{k). By the above observations on periodicity, and the fact that the path moves diagonally on 
Os, we deduce that 



efe+i 



+ efc. 



_ 2. 
The proof is concluded using induction and obvious monotonicity properties of the x-path. D 



Among many unresolved questions, we do not know whether an analogue of Proposition 16.5 
holds when the x-path is defined similarly in terms of the 1 Or 3 CA A rather than ^. 



7 Additive dynamics from random seeds 

Our goal in this section is to transfer the non-percolation results for infinite random initial 
configurations to random seeds. The proofs exploit an intriguing interplay between randomness 
and periodicity in the configuration started from a random seed. 

Lemma 7.1 (Random edge- intervals). Assume Aq is a uniformly random binary seed on [0, L]. 
For a fixed t, the state on the interval [t, t + L] 'x {t} is uniformly random. 

Proof. This is an immediate consequence of Lemma I3.2( i) and Corollary 14.31 D 
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Figure 7.1: Illustrations of the set of all points reached by paths starting within an initial random 
seed. Top: empty diagonal paths; bottom: wide paths. 

Lemma 7.2 (Edge-periodicity). For any Aq which is on [L + l,oo), and any k > 1, the 
sequence of edge configurations (A(i, t) : i = t + L — k + l,...,t + L) is periodic in t, with period 
at most 2k. 



Proof. This follows from Lemma l3.2r iv) and Lemma |4.1[ 



D 



Our first result establishes that, in subcritical cases, paths from the initial state do not reach 
far into the forward cone of [0, L] x {0}. This is illustrated in Figure [TTTl in which L = 25 and 
all points on paths from Z x {0} are again depicted in blue (only one layer of points outside the 
forward cone is colored blue, as all such are trivially reachable from Z x {0}). 

Proposition 7.3 (Percolation into the cone). Suppose Aq is a uniformly random binary seed on 
[0, L]. The probability that there is an empty diagonal path from Z x {0} to the forward cone of 
[0, L] X { [C log L\ } goes to as L —^ oo. The same is true for wide paths. Here C is an absolute 
constant. 



Proof. Let k he a positive integer to be chosen later satisfying 2k + 1 < L. Call a space-time 
point {x,t) bad if there exists an empty diagonal path from Z x {t — fc} to (x,f). If the state 
on the interval I{x, t) := [x — k, x + k] x {t — k} is uniformly random, then Theorem 11.51 implies 
that P((x,t) is bad) < exp(— cA;) for an absolute constant c > 0. 
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Figure 7.2: An illustration of the proof of Proposition 17.31 The outline of the forward cone of 
[0, L] X {0} is shown by the solid line (here L = 16). Points in the set S are shown as black 
discs. Any path from the top row to the region below S must pass through S. For each point 
in S, the top row of the associated triangle (of size A; = 3) has a uniformly random state (three 
such triangles are shaded). 

We define an infinite set of points S via Figure 17.21 This set has the following properties: 
(i) any path from Z x {0} to the forward cone of [0, L] x {2k} must pass through a point in S; 
and (ii) for every (x, t) € S, the state on the interval I{x, t) defined above is uniformly random, 
either trivially or by Lemma l7.1[ 

We wish to bound the probability that S contains a bad point by a union bound. The set S 
is infinite, but Lemma 17.21 implies that the states of the relevant intervals I{x,t) for (x,t) in the 
diagonal "arms" of S repeat with period at most 2(2A; + 2). Thus, besides the at most L points 
on the top section of S, there are only 8(2/c + 2) distinct cases to consider. Hence 



P{S contains a bad point) < {L + 16k + 16) e 



-ck 



The proof is completed by taking k = [C log L\ for a suitably large C (the argument for wide 
paths is identical). D 

Similarly, we next show that to each void of A* there corresponds a periodic strip that blocks 
diagonal and wide paths. Fix a void V of A', and an integer L > 0. Define the perturbed 
void V to be the triangular region 

WL{v) = vr^{v + {L,{))). 

See Figure [73] for an example. Note that Wiiy) = unless the width of V exceeds L. Further, 
fix an integer m > 1, assume the top of WLiV) = [a, 6] x {t}, and define the following interval 
above Wl{V): 

J^„(y) = [a- 2^,5 + 2^x^-2'". 

(We set JL,m{V) = when Wl{V) = 0.) 
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Figure 7.3: Perturbed void (dark blue), with L = 5, of the void with top interval [1, 15] x {16}. 
The perturbed void is filled with Os for any seed included in the interval of six red points. The 
forward cone of this interval is outlined. 
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Figure 7.4: Illustration of Theorem I7.4l fii) with the seed 110100111 on [0,8] and m = 2 and 
same two voids as in Figure [331 the repeating string is ^4 = 010110101111. 

Lemma 7.4 (Periodic and random intervals above voids). Suppose the initial configuration Aq 
vanishes outside [0,L]. Let m he a nonnegative integer. Let V be a void of X* of width at least 
L and at least 2™. 

(i) A vanishes on Wl{V). 
(ii) There exists a string A of length 3-2™', depending on m and Aq hut not on V, such that 

the configuration of X on JL,miy) is a suhword of A°^ . 
(iii) Now suppose that 2"^^^ < L and that Aq is uniformly random on [0, L]. Then every interval 

of length 2'" in JL,m(y) has uniformly random state. 



Proof. Claim (i) is a simple consequence of Lemma 14. 1^ (ii) follows from Proposition 13.61 and 
Lemma 14. H and (iii) from Proposition 13.61 and Corollary 14.31 D 
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Proposition 7.5 (Percolation into voids). Assume the conditions in Lemma \7.4^ iii). Let 
cross(m) be the event that there exists a void V for which there is either an empty diagonal or 
a wide path from JL,miV) to Wl{V). Then 

P(cross(m)) < exp(-c2'") 

for some universal constant c. 

Proof. Using Theorem 11.51 and Lemma 17.41 (ii) and (iii) , this follows by a similar argument to 
the proof of Theorem 17.31 The key point is that by Lemma I7.4l fii) , only 3 ■ 2"* distinct cases 
need to be considered in the union bound. D 

8 Replication and ethers in web cellular automata 

We can now prove Theorems 11.11 and 11.41 from the introduction. 

Proof of Theorem \1.4\ This is an immediate consequence of Proposition [73] and Lemma [2. 11 D 

Proof of Theorem \l.l\ Consider a uniformly random binary seed on [0,L]. In the context of 
Proposition 17.51 let M be the smallest m with 2*"^^ < L, for which cross(?ri) does not occur. 
If such an m does not exist, let M = oo. By Proposition 17. 5| M is tight as L — ?> oo. 

If M = oo, take Rl = oo. Assume now that M < oo. Then there exists a string A' of 
Os and 2s of length of 3 • 2™ so that the top row of Wl is a segment of {A')°° for every void. 
This holds because the first level configuration A is periodic with the required period on a strip 
above Wl, by Lemma 17.41 while the absence of relevant paths makes the top row also periodic 
by Lemma l2.1l Moreover, by the same results, the periodic pattern is the same for all voids. 

Consider the CA ^ started with a periodic configuration i?°°, for some string B of length 
a. The evolution is periodic in time after some initial burn-in period Tg. Let burnin((T) = 
maxsTe. Our random distance Rl is burnin(3 • 2 ) plus a universal additive constant. The 
proof is finished by Lemma 12.31 D 

As remarked earlier, Theorem 11.11 implies that the union of the regions that are filled by a 
translate of the ether has density 1 within the forward cone of the seed. Therefore, on the event 
that Rl < oo the set of non-zero points has a rational density within the same forward cone. 
We do not know whether the same holds for every seed. 

For an arbitrary web CA satisfying the conditions of Theorem ll.il limsup^_^oQ P(i?i > r) 
decays at least as fast as a power law in r. This is easily seen from the above proof, using the fact 
that burnin((T) < 2'^ . In cases when the dynamics restricted to Os and 2s is additive, including 
Web-Xor, Modified Web-Xor, and Piggyback, one can easily show that the decay is exponential. 
Identical remarks apply to the temporal period of the ether 7]l- 
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9 Bounds on ether probabilities 

In this section we prove Theorem 11.21 and explain how expUcit lower bounds on ether probabil- 
ities are proved. We also indicate how some ethers can be ruled out for certain rules. 

For ?Ti > 0, we call the string A in Lemma l7.4( ii) the level-2™ link of the seed Aq. (Note 
that the choice of A is unique up to periodic shifts.) Fix an integer fc > 1 and a binary string 
A. Consider 1 Or 3 with initial periodic configuration Aq = j4°°. If there is no empty diagonal 
(resp. wide, 6-iree) path from Z x {0} to Z x {/c — 1} in the resulting configuration A, then we 
say that that ^ is a blocker to depth k for diagonal (resp. wide, 0-free) paths. If Xk-i ^ 0, 
then we say that that A is non-degenerate to depth k. 

Fix an 771 > 1, and let Aq be a uniformly random seed on [0, L], with L > 2™-+^. For diagonal 
and wide paths. Proposition 17.51 implies that 

P(the level-2™ link of Aq is a blocker to depth 2™) > 1 - exp(-c2"), 

for some universal constant c. 

Further, consider a web CA ^t with an ether rj € {0, 2} . The signature of r; is a string 
B such that, for some t, r]{-,t) equals a (spatial) translation of B°°, and is the first in the 
lexicographic order among shortest such strings. Observe that two ethers are equivalent if and 
only if they have the same signature. We say that a binary string A produces r] with signature 
B if the initial state .^o = ^°° makes ^t equal to a translation of B°° at some time t. 

Lemma 9.1 (Blockers). Let S,t be a diagonal- compliant (respectively: wide- compliant, or 6- 
free- compliant) web CA. Further, let A be a string that is a blocker to depth 2™ for diagonal 
(respectively: wide, or 9-free) paths and produces an ether rj. If a seed S,o results in the level-2"^ 
link A, then ^ is a replicator with ether rj. If, in addition, the CA ^j has no spontaneous birth, 
then r] = 0. 

Proof. The first claim follows by the same arguments as in the proof of Theorem 11.11 The last 
claim follows by Lemma |2. 11 D 

Denote by Seedsj^^ the set of binary seeds that vanish outside [a, 6], and by Seeds = 
Ua<bSeeds[a ;,] the set of all binary seeds. Let g : {0, 1}^ — )• {0, 1}^ be the map determined by 
one step of the 1 Or 3 rule (i.e., the map Aq i-> Ai). It is well-known (and easy to prove) that, 
for a < b, the map g is injective from Seedsr^u to Seedsr^-i^fc+ii, and therefore the restriction 
fflseeds IS injective. We say that a binary seed Aq has a predecessor if it is in the image of 
ff I Seeds- More generally, Aq has k predecessors, for k > 1, if it is in the image of the A;-th 
iteration (5 1 seeds) ; in that case, (5 1 seeds )^ (Aq) contains a unique seed called the the k-th 
predecessor of Aq- We denote by Pred^ the set of all seeds that have k predecessors. The 
following lemma follows immediately from the properties of g' | seeds- 

Lemma 9.2 (Predecessors of random seeds). Assume Aq is a uniformly random binary seed on 
[0, L] and that 1 <k < L/2. Then P(Predfc) = 1/4 . Moreover, conditioned on Pred^, the k-th 
predecessor of Aq is a uniform binary seed on [k, L — k]. 
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Lemma 9.3 (Predecessors and links). For m >0, a seed Aq has 2™" predecessors if and only if 
its level-2'^ link is 0. 

Proof. If Aq has k predecessors, then its level-2™ link is by Lemma l3.2( ii) and Lemma l4.ll 

Conversely, assume that Aq G Seedsro^^] is given by a string S of length L + 1. If n is 
large enough so that 2*^ > 2L and 2^ > 2"^"*"^, the configuration of A2'i on [0,L] is S, again by 
Lemma [3.2( ii) and additivity. If the level-2™ link is 0, the configuration on [0, L] at time 2" — 2"^ 
provides a seed Ag, such that g'^'^iX'o) = Aq. D 



Lemma 9.4 (Ether probabilities). Suppose some seed Sq £ Seedsro^^-i] 'i'^ o, replicator with 
some ether rj, and is such that for some m > 1, the level-2"^ link is a blocker to depth 2'". Let 
(,0 be a uniformly random binary seed on [0, L] . Then 



liminf P(^o is a replicator with ether rf) > 2' 



-s-2^+1 



Proof. Assume a seed Si with support in [s, oo) has 2™ predecessors. Form a seed S by adding 
the configurations of So and Si. Then, by Lemma 19.31 and Lemma l4.lt S has the same level-2'" 
link as Sq, and therefore, by Lemma |9.H is a replicator with the same ether r/. In the rest of 
the proof we apply this fact to random seeds. 

Suppose now ^o is a uniformly random seed in Seedsrg^n) with L large enough so that 
L — s > 2™+-'^. Let ^Q (resp. ^g) be the random seed that agrees with ^q on [0, s — 1] (resp. [s, L]) 
and vanishes elsewhere. Then, 

P(5' is a replicator with ether r]) >P{(,q = So, and ^q has 2"* predecessors) 

= p(^^ = So) • P(Co has 2™ predecessors) 

= 2-^-4-2'", 

where the last equality follows from Lemma l9.2i D 

Proof of Theorem \1.^ This is immediate from Lemmas 19.11 and 19.41 and the proof of Theo- 
rem [TTTl n 



Recall that Extended 1 Or 3 is not diagonal- or wide-compliant. However, it is 4- free com- 
pliant, and this allows us to prove the following lower bounds. 

Theorem 9.5 (Replication and ether probabilities for Extended 1 Or 3). Let ^ be the Extended 
1 Or 3 web CA, started from a uniformly random binary seed on [0, L]. Then 

liminf P(^ is a replicator) > 0.826. 

L—>-OD 

Moreover, lower bounds on lim inf 2,_^oo P (C is a replicator with ether ij) for certain ethers rj are 
as in Table f^T/l 
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ether 
signature 


temporal 
period 


spatial 
period 


density 
of 2s 


lower 
bound 





1 


1 





0.6061 


02 


1 


2 


1/2 


0.0471 


0002 


2 


4 


1/2 


0.0333 


[7]2 


4 


8 


3/8 


0.0664 


5 202 


4 


8 


3/8 


0.0189 


15 2 


8 


16 


5/16 


0.0193 


[13]202 


8 


16 


11/32 


0.0079 


[11] 20002 


8 


16 


5/16 


0.0024 


[9] 2000202 


8 


16 


3/8 


0.0085 


[9] 2020202 


8 


16 


3/8 


0.0006 


[7] 200020202 


8 


16 


13/32 


0.0045 


[7] 202020202 


8 


16 


7/16 


0.0105 


[5] 2 [5] 20202 


8 


16 


7/16 


0.0006 



Table 9.1: Some non-equivalent ethers that provably emerge for Extended 1 or 3 from a long 
random seed with positive asymptotic probability. Each ether is generated from the initial 
condition obtained by repeating its signature indefinitely. Here [k] stands for an interval of k 
Os. The last column is a rigorous lower bound for the lim inf of the probability in Theorem 11.21 
The lower bounds sum to just over 0.826. 

For an ether rj, its reflection around the time axis is denoted by f/. Then r/ is symmetric 
if fj is equivalent to r/. Assume that B is the signature of ij and B its reflection. A sufficient 
condition for symmetry of r] is that the reflection i? is a periodic shift of B. However, this 
is not a necessary condition: the ether with signature B = [7] 200020202 is symmetric as the 
fourth iteration of the 1 Or 3 rule applied on B°° yields a translation of B°°. Thus the only 
non-symmetric ether in Table [9T] is the one with signature [9]2000202. In non-symmetric cases, 
our tables combine the frequencies of an ether and its reflection. 



Proof of Theorem \9.5[ Throughout the proof, fix a positive integer m > 1 and assume that 
L > 3 • 2™. Using the same notation as in Proposition 13. 6( assume that for some a > 0, the 
configuration of A* at some time ton/= [a-|-L,a-|-L — 1-I-3- 2™] x {t} is exactly the string 
^0 = 1 I I 1 1^^ I I . For ease of reference, we will assume A^ is positioned on Z so that ^o 
is the configuration in [0, 3 • 2" — 1]. 

Our main tool is the map <1> : Zg ^ — )• Zg''^™ that takes as argument an initial binary seed Aq 
supported on [0, L] and outputs the configuration of A on /. This is a linear map that assigns 
to every seed with support in [0, L] its level-2'" link. The matrix of $ (in the standard basis) 
has row i given by the segment [i, L + i] of A^^, i = 0, . . . , 3 • 2™ — 1. (All matrix and vector 
coordinate indices start at 0.) It is easy to see that the matrix has rank 2'""'"^, and therefore its 
image has cardinality 2^'" . The kernel of ^* has basis vectors y^, /c = 0, . . . , 2'" — 1, given by 
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y^ = l[i mod 2™ = A;], i = 0, . . . , 3 • 2"^ — 1. Therefore, the image of 'i' is the set 

$(Z^+i) = {be Zf^" : bi + b2m+i + 62^+i+i = 0, Vi = 0, . . . , 2"" - l}. 

A vector in Zg'^ is naturahy identified with a binary string of length 3 • 2*" and we will do so 
for the rest of the proof. Let A^^^ be the number of strings in ^{1^2 ) that are non-degenerate 
to depth 2™. Further, let A*";, the number of strings in ^{1^2^ ) that are non-degenerate and 
blockers, for 4-free paths, to the same depth 2™. 

Now suppose that ^o be a uniform random binary seed on [0, L] and let pL be the probability 
that ^ is a replicator. We claim that 

(9.1) liminfpi > — -. 

Recall that Predi is the event that ^o has a predecessor; by Lemma [9.31 Predf is exactly the 
event that $(.^o) is non-degenerate to depth 2™. Furthermore, conditioned on Predi, the first 
predecessor of ^o is a uniformly random binary seed on [1,L — 1]. Therefore, 

PL > P(^(^o) is a blocker to depth 2"* | Predf ) P(Predf ) 
/q 2\ + P(^ is a replicator | Predi) P(Predi) 

AT, 3 , 1 

Now (j9.ip follows by taking liminf as L — > oo of the first and last expressions of ()9.2p . The 
particular bound was obtained by a computer for m = 4: all 2^^ vectors in the range of ^ were 
checked for blocking and non-degeneracy, and the resulting tallies were A„ = 3,221,225,472 
and Nb = 2, 663, 229, 504. This concludes the proof for replication probability. 

The proof for a lower bound for a particular ether ij is identical except in the definition of 
Nh, which is now the number of strings in ^{'Z2~^ ) that are blockers and non-degenerate to the 
level 2™, and produce rj. For example, the result for the zero ether was A^f, = 1, 952, 489, 232. D 

Table 19.11 suggests that spatial and temporal periods of Extended 1 Or 3 ethers are powers 
of 2, and that the ether (2)°° never appears. This is addressed in our next two results. 

Lemma 9.6 (Periodic configurations). Assume that Aq is a spatially periodic configuration 
whose period a divides 3 • 2", and that Xt = Aq for some t. Then a divides 2". Moreover, for 
cr > 1 the temporal period equals a/2. 

Proof. By Lemma [3. II and Lemma [4. 11 we may assume n = 1, and then we check that any Aq of 
period 3 leads to a constant configuration in a single time step. The last assertion follows from 
Lemma l3. II and the following two easily checked facts: (1) if Aq is periodic with period at most 
2, then Aq = Ai; and (2) if Aq is periodic with period exactly 4, then Aq / Ai. D 

Proposition 9.7 (Possible ethers). Assume S,t is the Extended 1 or 3 CA. Suppose that that S,o 
is a replicator with ether rj, and that its level-2^ link is a blocker to depth 2"* for A-free paths. 
Then rj has spatial period that is a power of 2. Also, the signature of rj is either or it is of 
the form [ai]2[a2]2 . . . [ak\2, where k > 1 and each [a,,] is a string of Os of odd length ai. In 
particular, rj ^ 2. 
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ether 
signature 


temporal 
period 


spatial 
period 


density 
of 2s 


lower 
bound 





1 


1 





0.5 


2 


1 


1 


1 


0.0398 


02 


1 


2 


1/2 


0.0142 


0002 


2 


4 


1/2 


0.0258 


[7 2 


4 


8 


3/8 


0.0099 


[4 2022 


4 


8 


1/2 


0.0303 


00020222 


4 


8 


1/2 


0.0209 


00022222 


4 


8 


5/8 


0.1297 


0002000200022222 


8 


16 


11/16 


0.0362 


0002000200202002 


8 


16 


9/16 


0.0216 



Table 9.2: Ten ethers that emerge from long random seeds for Piggyback with positive asymptotic 
probability. The conventions of Table 19.11 apply. 

Proof. The first claim follows from the previous lemma and Theorem 17.41 so we proceed to 
prove the second claim. If At = 0, but Xt-i ^ 0, then there are, up to translation, exactly two 
possibilities for Ai_i and Xt-2- 



111100111100 ■•• 
011011011011 ••• 



•••010001010001 
•••011011011011 



Assume that the seed ^o is such that S{^o) has exactly k predecessors. (Recall that 5{a) := 
l[a = 1].) Then, for any n, the state of 6{C) on [C - 1, 2" - C + 1] x {2'^ - A; - 2, 2" - fc - 1} is a 
segment of one of the two configurations above. (Here, C is a constant that depends only on L.) 
By considering 4-free paths, we see that the left configuration implies ^2"-fc = on [C, 2" — C], 
while for right one implies that one ^2"-A; vanishes outside 6Z n [C, 2" — C] (after a suitable 
translation). As 2s evolve according to the 1 Or 3 rule in the absence of Is, the positions of 2s 
started from a subset of 2Z are a subset of 2Z at all even times (by Lemma l3.ip . The claimed 
form of the signature follows. D 



Lower bounds for ether probabilities can also be obtained for Piggyback, with the same proof 
as for Theorem 19^ 



Theorem 9.8 (Ether probabilities for Piggyback). Let ^t be the Piggyback web CA, started from 
a uniformly random seed ofOs and Is on [0,L]. Lower bounds on 



liminf P(^ is a replicator with ether rf) 



are as in TablelKE 



The computer search with m = 4 yielded 117 different ethers for Piggyback with provably 
positive asymptotic probability, with their combined probabilities at least 0.914. The ethers 
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listed in Table [9^2] are the ten with largest lower bounds. The only non-symmetric ether among 
these ten has signature [4]2022. (The initial state (00020222)°° generates its translated reflection 
in 2 steps.) We do not know whether the asymptotic probability for the zero ether is exactly 
1/2. 

Open problems 

As the earlier discussions indicate, this topic offers a rich supply of open questions. We highlight 
a small selection. 

(i) In the 1 Or 3 cellular automaton started from a uniformly random binary string on the 
half-line [0,oo), what is the growth rate of the maximum integer rt for which there is an 
empty path from (— oo,0) x {0} to (rt,t)? Is it the case that rt/t —^ 1/4 as t — )• oo? 

(ii) What can be said about percolation in the space-time configuration of other one- 
dimensional cellular automata started in an invariant measure? For example, the uniformly 
random binary string on Z is invariant for permutative rules (see jBLj for a definition), 
including Rule 30. Do there exist infinite diagonal, wide or empty paths? 

(iii) Are there infinitely many different ethers for replicators in the Piggyback cellular automa- 
ton? Is there an algorithm that decides whether a given ether occurs in some replicator? 

(iv) For two-dimensional Box 13 solidification CA (see Section [2] and |GG2| ) started from a 
uniform random seed in [0, L]^, does the final configuration have rational density with 
probability converging to 1 as L — )■ oo? 



Acknowledgements 

JG gratefully acknowledges the hospitality of the Theory Group at Microsoft Research, where 
most of this work was completed. He was partially supported by NSF grant DMS-0204376 and 
the Republic of Slovenia's Ministry of Science program P 1-285. 



References 

[BBS] P. N. Balister, B. Bollobas, A. M. Stacey, Dependent percolation in two dimensions, 
Probability Theory and Related Fields 117 (2000), 495-513. 

[BDR] C. D. Brummitt, H. Delventhal, M. Retzlaff, Packard snowflakes on the von Neumann 
neighborhood, Journal of Cellular Automata 3 (2008), 57-80. 

[BS] R. Basu, A. Sly, Dependent percolation in two dimensions, larXiv: 1204.29311 (2012). 

[BL] M. Boyle, B. Lee, Jointly periodic points in cellular automata: computer exploration and 
conjectures, Experimental Mathematics 16 (2007), 293-302. 



REFERENCES 43 

[BBDM] J. Balogh, B. Bollobas, H. Duminil-Copin, R. Morris, The sharp threshold for bootstrap 
percolation in all dimensions, Transactions of the American Mathematical Society 364 
(2012), 2667-2701. 

[Coo] M. Cook, Universality in elementary cellular automata, Complex Systems 15 (2004), 1-40. 

[EN] K. Eloranta, E. Nummelin, The kink of cellular automaton Rule 18 performs a random 
walk. Journal of Statistical Physics 69 (1992), 1131-1136. 

[EGG] R. Fisch, J. Gravner, D. Griffeath, Metastability in the Greenberg-Hastings Model, Annals 
of Applied Probability 3 (1993), 935-967. 

[Gael] P. Gacs, Reliable cellular automata with self- organization. Journal of Statistical 
Physics 103 (2001), 45-267. 

[Gac2] P. Gacs Compatible sequences and a slow Winkler percolation, Combinatorics, Probabil- 
ity and Computing 13 (2004), 815-856. 

[GGl] J. Gravner, D. Griffeath, Modeling snow crystal growth I: rigorous results for Packard's 
digital snowflakes. Experimental Mathematics 15 (2006), 421-444. 

[GG2] J. Gravner, D. Griffeath, Asymptotic densities for Packard Box rules, Nonlinearity 22 
(2009), 1817-1846. 

[GG3] J. Gravner, D. Griffeath, The one- dimensional Exactly 1 cellular automaton: replication, 
periodicity, and chaos from finite seeds, Journal of Statistical Physics 142 (2011), 168-200. 

[GG4] J. Gravner, D. Griffeath, Robust periodic solutions and evolution from seeds in one- 
dimensional edge cellular automata. Theoretical Computer Science 466 (2012), 64-86. 

[GGP] J. Gravner, G. Gliner, M. Pelfrey, Replication in one- dimensional cellular automata, 
Physica D 240 (2011), 1460-1474. 

[Gra] L. F. Gray, A Reader's Guide to Cacs's 'Positive Rates" Paper, Journal of Statistical 
Physics 103 (2001), 1-44. 

[Hoi] A. E. Holroyd, Sharp metastability threshold for two-dimensional bootstrap percolation. 
Probability Theory and Related Fields, 125 (2003), 195-224. 

[Lin] D. A. Lind, Applications of ergodic theory and sofic systems to cellular automata , Physica 
D: Nonlinear Phenomena 10 (1984), 36-44. 

[MOW] O. Martin, A. M. Odlyzko, S. Wolfram, Algebraic properties of cellular automata. Com- 
munications in Mathematical Physics 93 (1984), 219-258. 

[Pil] S. Y. Pilyugin, Shadowing in Dynamical Systems, Lecture Notes in Mathematics Vol. 1706, 
Springer- Verlag (1999). 

[Wil] S. J. Willson, Cellular automata can generate fractals, Discrete Applied Mathematics 8 
(1984), 91-99. 



REFERENCES 44 



[Win] P. Winkler, Dependent percolation and colliding random walks, Random Structures & 
Algorithms, 16 (2000), 58-84. 

[WoU] S. Wolfram, Statistical Mechanics of Cellular Automata, Reviews of Modern Physics 55 
(1983) 601-644. 

[Wol2] S. Wolfram, Random sequence generation by cellular automata. Advances in Applied 
Mathematics 7 (1986), 123-169. 



